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A search for the relic neutrinos from all past core-collapse supernovae was conducted using 1496 
days of data from the Super-Kamiokande detector. This analysis looked for electron-type anti- 
neutrinos that had produced a positron with an energy greater than 18 MeV. In the absence of a 
signal, 90% C.L. upper limits on the total flux were set for several theoretical models; these limits 
ranged from 20 to 130 i^e cm~'^ s~^. Additionally, an upper bound of 1.2 Ue cm~^ s~^ was set for 
the supernova relic neutrino flux in the energy region > 19.3 MeV. 

PACS numbers: 95.85.Ry,97.60.Bw,14.60.Lm,96.40.Tv 



During a core-collapse supernova, approximately 10^"^ ergs of energy are released, about 99% of which are in the 
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form of neutrinos. To date, the only time that a burst 
of such neutrinos has been detected was in the case of 
SN1987A |. However, it is generally believed that 
core-collapse supernovae have occurred throughout the 
universe since the formation of stars. Thus, there should 
exist a diffuse background of neutrinos originating from 
all the supernovae that have ever occurred. Detection 
of these supernova relic neutrinos (SRN) would offer in- 
sight about the history of star formation and supernovae 
explosions in the universe. 

All types of neutrinos and anti-neutrinos are emitted 
from a core-collapse supernova, but not all are equally 
detectable. The Df. is most likely to be detected by Super- 
Kamiokande (SK). It interacts primarily through inverse 
(3 decay {ve P ^ n , Ee = - l.i MeV) with a 
cross section that is two orders of magnitude greater than 
that of neutrino-electron elastic scattering. All further 
discussion herein of the SRN refers only to the ■ 

Several methods have been used to model the SRN flux 
and spectrum j|, |[ ||, 0, H) , with flux predictions rang- 
ing from 2-54 cm~^ s . In this paper, SK search re- 
sults are compared to SRN predictions based on a galaxy 
evolution model Q , a cosmic gas infall model [|j , cosmic 
chemical evolution studies ||^ , and observations of heavy 
metal abundances ||^. A model that assumes a constant 
supernova rate Q was also considered; this model was 
used to set the previous SRN flux limit at Kamiokande- 

n [|. 

It has been shown that neutrinos undergo flavor os- 
cillation 0, |ll|, H). Therefore, an SRN spectrum that 
includes the effects of neutrino mixing was also consid- 
ered All six species of neutrinos are emitted during 
a core-collapse supernova. However, the and the Ur 
decouple from the ncutrinosphere earlier than the v^, re- 
sulting in a higher temperature for these flavors; thus, 
neutrino mixing would harden the energy spectrum. 
The cross section for inverse /3 decay increases as the 
square of the Ve energy, so oscillation would enhance the 
SRN signal. In this paper, only a large mixing angle 
MSW solution (LMA) was used to distort the SRN spec- 
trum - the LMA solution is also favored by the available 
solar neutrino data . 

This paper presents the results of a search for SRNs 
in the Super- Kamiokande detector. SK is a water 
Cherenkov detector, with a fiducial mass of 22.5 kton, 
located in the Kamioka Mine in Gifu, Japan. Descrip- 
tions of the detector can be found elsewhere |l5|. The 
data reported here were collected between May 31, 1996, 
and July 15, 2001, yielding a total SRN search livctimc of 
1496 days. Backgrounds to the SRN signal are solar neu- 
trinos, atmospheric neutrinos, and muon- induced spalla- 
tion products. Background reduction takes place in the 
following steps: spallation cut, sub-event cut, Cherenkov 
angle cut, and solar direction cut. 

Spallation is the most serious background, and the 
ability to remove it determines the lower threshold of 
the SRN search. Spallation products are also relevant to 
solar neutrino studies, and so a likelihood function had 



been developed that uses information about the muons 
preceding the possible spallation event [ll| . To permit a 
low analysis threshold for the SRN search, a tighter spal- 
lation cut was implemented: in addition to the likelihood 
function cut, all events that occur less than 0.15 seconds 
after a cosmic ray muon are rejected. The spallation cut 
is applied to all events reconstructed with E < 34 MeV, 
and introduces a deadtime of 36%. No discernible spal- 
lation events with energies above 18 MeV remain in the 
data after this cut is applied and so 18 MeV was set as 
the lower analysis threshold. 

The sub-event cut removes muons with kinetic energy 
(T) in the range of 50 - 140 MeV. Cosmic ray muons, 
which have much higher energies, originate outside of SK 
and are removed because they produce a veto signal in 
the SK outer detector (OD). However, atmospheric z/^ 
can produce muons within the inner detector (ID) via 
charged current quasi-clastic scattering {v^ N — > /x N')] 
such muons will not have been tagged by the OD, but 
will be visible in the ID. Muons with low energies will 
stop in SK and produce a decay electron; often the muon 
and decay electron are found in the same event. When 
this happens, the decay electron is referred to as a "sub- 
event." After the vertex of each event was found and the 
flight time of the Cherenkov photons was subtracted, the 
event's 1.3 /xs time window was searched; if more than one 
timing peak was present, then the event was removed. 
The sub-event cut was tested on simulated muons in the 
relevant energy range, and it was shown to remove 34% 
of the muon background. 

The remaining muons are removed by the Cherenkov 
angle cut. This cut exploits the mass difference be- 
tween the muon and the positron, which results in a 
difference in their Cherenkov angles 9c- Positrons with 
E > 18 MeV have Qc ~ 42 degrees. Muons with 
T < 140 MeV have Oq < 34 degrees; thus, all particles 
with 9c < 37 degrees were removed from the data. The 
efficiency of this selection criterion for retaining signal is 
98%, as determined by applying it to simulated positron 
events. Using simulated muon events, it was shown that 
applying the Cherenkov angle cut and the sub-event cut 
together results in the rejection of 96% of the muon back- 
ground. Furthermore, it was found that the full reduction 
removes > 99% of the muons. The Cherenkov angle cut 
was also used to remove events with 9c > 50 degrees; this 
eliminated events without clear Cherenkov rings, such as 
multiple 7 rays emitted during a nuclear de-excitation, 
from the data sample. 

Finally, a cut on the direction of the event is made to 
remove contamination from solar neutrinos. Events with 
E < 34 MeV were removed if the reconstructed event 
direction was within thirty degrees of the vector from 
the Sun to the Earth at the time of the event. 

By simulating positrons created from SRN, it was 
found that the efficiency of the full data reduction is 
47±0.4% for < 34 MeV, and 79±0.5% for S > 34 MeV. 
Figure [l] plots the energy spectrum after each cut. 

After applying the selection criteria, two irreducible 
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FIG. 1: Energy spectrum at each reduction step. In the final 
data set, the spaUation cut and solar direction cut are only 
applied in the first four bins. The numbered lines represent 
the corresponding theoretical predictions from Table M. 



FIG. 2: Energy spectrum of SRN candidates. The dotted and 
dash-dot histograms are the fitted backgrounds from invisible 
muons and atmospheric . The solid histogram is the sum of 
these two backgrounds. The dashed line shows the sum of the 
total background and the 90% upper limit of the SRN signal. 



backgrounds remain. The first is atmospheric and 
j/g events. The second comes from atmospheric that 
interact to form a muon with T < 50 MeV. The en- 
ergy of these muons is below the threshold for emitting 
Cherenkov photons, so they are said to be invisible. De- 
cay electrons from visible muons can be eliminated; how- 
ever, when an invisible muon decays there is no way to 
tag the resulting electron as a background event. 

The energy spectra of these backgrounds have shapes 
that are very different from each other and from the SRN 
signal shape. In the region where SRN events are ex- 
pected (18 ~ 34 MeV), the dominant background is de- 
cay electrons from invisible muons, which have energies 
that are distributed according to the Michel spectrum. 
However, at higher energies, atmospheric Ve events dis- 
tort the Michel spectrum. To evaluate the distortion, it is 
necessary to extend the upper analysis threshold to ener- 
gies where only atmospheric v^, events are present. Decay 
electrons have a maximum energy of 53 MeV, but may 
be detected up to ^ 65 MeV due to the energy resolution 
of SK. Beyond 65 MeV, only atmospheric v^^ are found, 
so the upper analysis threshold was set above 65 MeV 
and the data were analyzed with a three parameter fit. 

To determine the final shape of the backgrounds, 
100 years of simulated events were generated per back- 
ground. The initial shape of the decay electrons was 
determined by the Michel spectrum; the initial shape of 
the atmospheric events was obtained from previous 
works ||l^, The background simulations were sub- 

jected to the full reduction, and the shape of the resulting 
spectra were used to fit the data; each of the SRN mod- 



els was treated similarly. For the fitting, the data were 
divided into sixteen energy bins, each 4 MeV wide (see 
Figure and the following 'x^ function was minimized 
with respect to a, /3, and 7. 

2 _ [(q ■ M) + (/3 ■ Bi) + (7 ■ Ci) - Nif 

In this equation, the sum I is over all energy bins and 
iV; is the number of events in the bin. Ai, Bi, and C; 
represent, respectively, the fractions of the SRN, Michel, 
and atmospheric spectra that are in the l*^ bin. a, 13, 
and 7 are the fitting parameters for the number of SRN 
events, decay electrons, and atmospheric Ve events. 

The total number of events in the data sample is small, 
so the statistical error Gstat is the dominant term in the 
denominator. The systematic error Usys considers the 
effects that uncertainties in the spectrum shapes have 
on the SRN result. Such uncertainties originate from 
the reduction, the SK energy resolution, the theoretical 
atmospheric v,, spectrum, and other sources. For all bins, 
CTsj/s ~ 6%, which is always much smaller than astat- 

The efficiency-corrected event rate spectrum of SRN 
candidates and the results of the fit are shown in Figure]^. 
The best fits to 7 and (3 are indicated, respectively, by 
the dot-dashed and dotted lines. The solid line is the 
sum of these lines and represents the total background. 
For all six models, the best fit to a was zero and the 
minimum value was 8.1 for 13 degrees of freedom. A 
90% C.L. limit on a was set for each model; the dashed 
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TABLE I: The SRN search results are presented for six theoretical models. The first column describes the method used to 
calculate the SRN flux. The second column shows the efficiency-corrected limit on the SRN event rate at SK. The third column 
is the flux limit set by SK, which can be compared with the theoretical predictions that are shown in the fourth column. The 
flfth column shows the flux predictions above a threshold of Ei, > f 9.3 MeV . Note that the heavy metal abundance calculation 
only sets a theoretical upper bound on the SRN flux M. 



Theoretical model 


Event rate limit 


SRN flux limit 


Predicted flux Predicted flux 




(90% C.L.) 


(90% C.L.) 


(E, > 19.3 MeV) 



Galaxy evolution \i\ 
Cosmic gas infall [5 
Cosmic chemical evolution 
Heavy metal abundance 
Constant supernova rate 
Large mixing angle osc. 



< 3.2 events/year < 130 i^e cm~ s~ 

< 2.8 events/year < 32 Ue cm~^ s~^ 

< 3.3 events/year < 25 Ue cm~^ s~^ 

< 3.0 events/year < 29 Ue cm^^ s' 

< 3.4 events/year < 20 Ue cm~^ s" 

< 3.5 events/year < 31 cm~^ s~ 



8.3 Us cm 



52 i^e cm 



0.41 Ue cm 
0.20 Ve cm" 



0.43 Ve cm 



line represents the sum of the background and the upper 
bound on a for the galaxy evolution model. This line 
shows the type of distortion in the Michel spectrum that 
would be indicative of an SRN signal. 

Figure || shows that the expected backgrounds fit the 
data well. In this analysis, no flux normalization was cho- 
sen for the background rates; only the shapes were used. 
This is because there are large uncertainties (~ 30%) in 
the atmospheric neutrino fluxes at these very low ener- 
gies. As a consistency check, the fit results for the num- 
ber of background events were compared to the predic- 
tions, which were determined by applying the reduction 
cuts to 100 years of simulated background events and 
normalizing for livetime. For 1496 days of data, the ex- 
pected number of atmospheric Ve events is 75 ± 23, which 
is consistent with the best fit result of 88 ± 12 events. 
To determine the expected number of decay electrons 
from invisible muons, neutrino oscillation must be con- 
sidered JlOt ; given the low energy of the atmospheric 
that produce invisible muons, it is assumed that half of 
the i^fi have oscillated into i^r- With this assumption, the 
predicted number of decay electron events is 145 ± 43, 
which is consistent with the best fit result of 174 ± 16 
events. 

The limit on a can be used to derive a 90% C.L. limit 
on the SRN flux from each model. The number of SRN 
events is related to the total flux F by the following equa- 
tion: 



In this equation, Np is the number of free protons in SK 
(1.5 X 10^^), T is the detector livetime (1496 days), e{E) 
is the signal detection efficiency, (j{E) is the cross section 
for the inverse (3 decay (9.52 x lO"^"* Eo Po), and f{E) is 
the normalized SRN spectrum shape. The integral spans 
the energy range of the neutrinos that produce positrons 
in the observed region. 



Using the above values, the 90% C.L. SRN flux limit 
was calculated for each model. The results are in the 
third column of Table |, and can be compared with 
the predictions, which are in the fourth column. For 
the galaxy evolution model the cosmic gas infall 
model and the cosmic chemical evolution model 
the SK limits are larger than the predictions by a factor 
of three to six. In these models, the dominant contribu- 
tion to the SRN flux comes from supernovae in the early 
universe, so the neutrino energy is red-shifted below the 
18 MeV threshold. The heavy metal abundance model 
primarily considers supernovae at red-shifts z < 1, so SK 
is sensitive to more of the SRN flux. For this model, 
the flux limit is smaller than the calculated total flux . 
However, this prediction is only a theoretical upper limit, 
so these results can constrain this model but they cannot 
eliminate it. The LMA model [|| has a harder energy 
spectrum, and so SK is sensitive to a larger fraction of 
the SRN flux. The increased sensitivity is offset by the 
fact that this hardened spectrum also results in a larger 
limit for a\ thus, the SK flux limit is still nearly a factor 
of three larger than the prediction. 

The total SRN flux predicted by the constant model 
scales with the rate of core-collapse supernovae, and so 
the SRN flux limit can be used to set a 90% C.L. up- 
per limit on the constant supernova rate. The SRN flux 
prediction quoted in this paper is based on a reasonable 
supernova rate of 1.6 x 10'^ SN year~^ Mpc"'^. The ob- 
served SRN flux limit (20 cm^^ s^^) corresponds to 
a supernova rate limit of 6.2 x 10^ SN year~^ Mpc~^. 
Thus, the constant model can be ruled out, as the limit 
on the supernova rate is too low to be consistent with the 
observed abundance of oxygen [|[ which is synthe- 
sized within the massive stars that become supernovae. 
At Kamiokande-II, a flux limit of 780 cm~^ s~^ was set 
with the assumption of a constant supernova model 
the SK limit is 39 times more stringent. 

The SRN limits vary greatly, based on the shape of 
the theoretical SRN spectrum at energies that are be- 
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low SK's SRN analysis threshold. To remove this strong 
model dependence, a limit was set for E^, > 19.3 MeV. In 
this region, all six models have similar energy spectrum 
shapes, and so an experimental limit that is insensitive 
to the choice of model can be obtained as follows: 



ng.sMoV 



f{E,)dE, 



Jo" f{E,)dE, 



(3) 



Flux limits in this energy region were the same for 
all models considered: 1.2 Vf, cm~^ s~^. Previously, the 
best limit on the SRN flux in this region was set using 
357 days of data at Kamiokande-II This limit was 
226 Ve cm~^ s~^; the current SK limit is two orders of 
magnitude lower. 

In summary, a search was conducted at Super- 



Kamiokande to detect the diffuse signal of Ve from all 
previous core-collapse supernovae. No appreciable signal 
was detected in 1496 days of SK data. Using various 
models, 90% C.L. limits were set on the total SRN flux. 
A limit of 1.2 cm~^ s"'^ was set for the SRN flux 
above a threshold of E^, > 19.3 MeV. These results are 
more than an order of magnitude better than previous 
limits; some theories regarding the supernova rate in the 
universe can be constrained or rejected by these limits. 
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